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Abstract We expand the equations used in MacDonald's 1964 theory and Fourier analyze the tidal 
variations of the height at one point on the Earth surface, or, alternatively, the tidal potential at such 
point. It is shown that no intrinsic law is relating the lag of the tide components to their frequencies. In 
other words, no simple rheology is intrinsically fixed by MacDonald's equations. The same is true of the 
modification proposed by Singer(1968). At variance with these two cases, the modification proposed by 
Williams and Efroimsky (2012) fix the standard Darwin rheology in which the lags are proportional to 
the frequencies and their model is, in this sense, equivalent to Mignard's 1979 formulation of Darwin's 
theory. 
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1 Introduction 

MacDonald's name is often used to designate tidal theories in which the tide lags are frequency inde- 
pendent (see e.g. Ferraz-Mello et al. 2008 footnotes 4 and 6). This is however wrong. Indeed MacDonald 
(1964) adopted a constant lag S and integrated the evolution equations keeping this value unchanged 
during the whole evolution. However, what is called lag in MacDonald's theory is not the same quan- 
tity as adopted in Kaula(1964) and/or in modern theories. Using Kaula's notations (see Williams and 
Efroimsky, 2012), frequency independent lags actually mean that all Si mpq are equal and independent 
of the frequencies ximpq- In Ferraz-Mello et al.'s notations, it means eq = sx = £2 = ■ • •. This does not 
happen in MacDonald's tide and thus it is not correct to associate Mac Donald's name to such tidal 
theories. 

We investigate in this note how lag the tide components in MacDonald's model. First we have to 
understand what is done in that paper. The key equations are MacDonald's eqns.(20— 21). Eqn. (20) 
is introduced as follows: 

Let us suppose that at the earth's surface the tide lags by an angle 8 , so that, at a time / 
(the moon's angular position, / , is used to keep time), the potential is not W(/) but rather 
W(f — 8). Thus, at a time / , the extra potential due to the tidal deformation of the surface of 
the earth is: 

W = k 2 W 2 (f - S) (MD. eqn. 20) 

since the rotation carries the bulge forward and the zenith distance becomes 5 
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(We omit the terms with subscripts 3, 4, . . . , n, because they are irrelevant in this discussion). 

This quotation is the point where ambiguity is introduced in MacDonald's paper because / is the 
true anomaly, which is not a uniform quantity, and which could not be used as equivalent to time. 

The only way to overcome this ambiguity is learning from MacDonald's equations what he actually 
did. This is done in Sections 2 and 3 of this note. The results of these two sections are founded on 
MacDonald's eqn. (21), which is given in his paper as follows: 



The external potential at a time / is 



GmR 5 

U = 3 3 ® k 2 P 2 (cos5) (MD. eqn. 21) 



r f r f- 



where r/ corresponds to the instantaneous distance to the moon. 

(again, we omit the terms with subscripts 3, 4, . . . , n; we also changed the notation used for the Legendre 
polynomials to the standard one). In this equation, r/s is the distance to the Moon when the true 
anomaly of the Moon is / — 5. 

In the several instances in which these equations are used, / is always the true anomaly. This can 
be clearly seen, for instance, in MacDonald's eqn. (50) giving the quantity C = rj^ 5 rj 4 . 



2 Tide heights in MacDonald's Earth 

We may solve the ambiguity of MacDonald's written text without introducing any interpretative hy- 
pothesis. MacDonald's eqn. (21) is the classical expression of the potential created by a prolate spheroid 
on an external point (the Moon) and we may solve the inverse problem: 

Which is the spheroid creating on the Moon the additional potential given by MacDonald's 
equation (21) ? 

Thereafter, for the sake of avoiding embedding the point that we want to elucidate into a huge 
amount of algebra, we consider the planar approximation in which the reference plane is the plane of 
the lunar motion, the origin axis passes through the Moon's orbit perigee and the Earth rotation axis 
is perpendicular to the Moon's orbital plane. None of these assumptions interferes with the nature of 
the problem being considered. 

The major principal axis of the spheroid may lie on the intersection with the reference (orbital) 
plane, of one cone whose axis is directed to the Moon and having half-opening 5. Mathematically we 
have two solutions, but one of them is on the wrong side (the tide is advancing the Moon's position 
instead of lagging it). The height over the sphere of the points of the spheroid creating the potential 
U are given by: 

5i? 4 

5r = — -J^fc 2 (2cos2<f + 1) (1) 
12r}_ s 

where *P is the angular distance from the generic points of the prolate spheroid to its major principal 
axis. If we consider this point in the equator of the body: & = <fi— (w + f — 5) or <P = (p — f + 5 (since we 
adopted ui = 0), where <p is the longitudfe of the point. The harmonic decomposition of the resulting 
tidal height at that point is 

5i? 4 / 

5r = j^k 2 (4 + 12 cos(2£ -2(f)- 26) + 12e cos(^ -6) + 18e cos(£ -2<j>-6)- 
24e cos(£ - 2(f> - 25) + 24e cos(3^ - 2<j> - 25) + 18e cos(3^ - 2<j> - 2,5) + 
18e 2 - 12e 2 cos(<5) + 9e 2 cos(20) - 18e 2 cos(2^ + 5) + 9e 2 cos(20 + 25) + 12e 2 cos(2£ -5) + 
6e 2 cos(2£ - 25) + 18e 2 cos(2^ -2<f>-5) + 6e 2 cos(2^ -2(f)- 25) - 54e 2 cos(2^ -2<f>- 35) + 
39e 2 cos(4<? - 2<j> - 25) + 54e 2 cos(4^ - 2<f> - 3<5) + 9e 2 cos(4<? - 2(f> - 45)) + 0(e 3 ). 

(2) 
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Since the Earth is rotating with angular velocity Q 1 we may write <j> = <f>o + fit. The simple 
inspection of the above equation shows that in the strict MacDonald model the Fourier decomposition 
of the height of the tide does not show any simple law relating the lags and the frequencies of the tide 
components. In particular, it is not true that in a MacDonald tide all tide components have the same 
lag. 

The conclusion is that, as MacDonald fixed the vertex of the Earth tide at a constant distance 
5 from the Moon, his model leads in a univocally determined way to a trigonometric series in which 
the components lag on a very particular way. Even the simple semi-diurnal component has parts 
with different lags: cos(2£ — 2(f) — 25), e 2 cos(2^ — 2(f) — 5), e 2 cos(2^ — 2(f) — 35) (which have the same 
frequency, but the lags 26, 6 and 36 respectively). The same behavior may be verified in higher orders 
of approximation. 



3 MacDonald's potential of the deformed Earth 

Another possibility is to consider MacDonald's Eqn. 21 as the additional potential at a point whose 
distance to the Earth is the distance of the Moon, r / , and whose angular distance to the vertex of the 
tidally deformed Earth is 6. This means that MacDonald indeed introduced a lag in the true anomaly 
and not in time. 

MacDonald's eqn. (21) can be easily converted into the equation for the potential in a generic point 
P(P,# 

GmR% , 

W^^K^) (3) 

P r f-8 

where 'P is the angular distance from the generic point to the vertex of the tidally deformed Earth. 
This is a classical equation. 

This equation for the potential is the only possible one which is consistent with MacDonald's Eqn. 
(20) and with the words highlighted after that equation, "the zenith distance becomes <5", in which 
the distance of the Moon to the vertex of the Earth's tidal prolate spheroid is a constant 6. 

In the following, for simplicity, we assume the point P (p,(f>) on the Earth's equator (which again 
is assumed to coincide with the plane of Moon's orbit). Then *P is again the difference between the 
longitude of the point, <j>, and the longitude of the delayed Moon, lo + f — 6. 

Therefore, 

U=^PP L k 2 P 2 (co S (<j)-f + 6)) (4) 

P r f-S 

where, for simplicity, we have again adopted oj = 0. This function can be expanded in Fourier series 
to give: 

GmR b / 

U = 16 3fl f fc 2 (4 + 12 cos(2l - 2(f) - 26) + 12ecos(^-(5) + 18e cos(£ - 2(f) - 5)- 
24e cos{£ -2(f)- 26) + 24e cos(3£ -2(f)- 26) + 18e cos(3£ -2(f)- 36) + 

18e 2 - 12e 2 cos(c5) + 9e 2 cos(20) - 18e 2 cos(2</> + 6)+ 9e 2 cos(2^ + 26) + 12e 2 cos(2^ -6)+ ( 5 ) 
6e 2 cos(2^ - 26) + 18e 2 cos(2£ -2(f) -6) + 6e 2 cos(2^ -2<j>- 25) - 54e 2 cos(2£ -2(f)- 35) + 
39e 2 cos(4^ -2(f)- 25) + 54e 2 cos(4£ -2(f)- 35) + 9e 2 cos(4£ -2(f)- 4(5)) + C(e 3 ). 

the periodic part of which is the same as found in Sec. 2. The Earth potential in the considered point 
is a Fourier oscillations packet due to the Moon raised tide. As discussed in the previous section, the 
potential due to the strict MacDonald model is formed by oscillations that lag in a way not correlated 
to their frequencies. 

The coincidence of this result with that of the previous section proves that, contrarily to what can 
be read in his written text, MacDonald has indeed introduced a lag in the true anomaly and not in 
time. 



4 



4 The Williams-Efroimsky modified theory 

Williams and Efroimsky (2012) consider that a delay cannot be introduced in the true anomaly, which 
is not uniform. In their paper, the delay is introduced in the mean anomaly, and MacDonald's eqn. 
(21) is rewritten in a form equivalent to 

what means that they are considering not the action of the Moon when its true anomaly was / — 5, 
but when its mean anomaly was I — 5. 

We may perform the same calculations as above, to obtain for one generic point P(p, <f>): 

U = Hr fc 2 (2 + 6 cos(2£ - 20 + 25) + 6e cos(^ + 6)-3e cos(£ -20 + 5) + 

d>p 6 a 6 \ 

21e cos(3^ - 20 + 3<S) + 3e 2 + 9e 2 cos(2£ + 25) - 15e 2 cos(2^ - 20 + 25) + ( 7 ) 
51e 2 cos(4£ - 20 + 4<5)) + 0(e 3 ) 

Notwithstanding the much more simple expression, a law relating lags and frequencies is not yet 
visible. For that sake, one more hypothesis is needed (which is the same operational hypothesis intro- 
duced by Mignard(1979) in his concise presentation of Darwin's theory): At a given time, the potential 
on one point on the Earth is equal to the potential which would be acting on it at a time r before the 
given time, would the Earth be inviscid. This means that we have to introduce not the longitude of 
the point but — Or, or — (S75/n). The above expansion then becomes 



U 



GmRa 



■k 2 (2 + 6 cos(2£ - 20 + 25(1 -f2/n)) + 6e cos(£ + 5) - 3e cos(^ - 20 + 5(1 - 2Q/n) + 



8p 3 a 3 

21e cos(3£ - 20 + 3(5(1 - 2J?/3n)) + 3e 2 + 9e 2 cos(2^ + 25)) - 15e 2 cos(2£ - 20 + 25(1 -Q/n)) + 
hie 2 cos(U - 20 + A5(l - 2/2/4n))) + <D(e 3 ) 

where Darwin's law of lags proportional to frequencies emerges making the new theory consistent with 
the standard Darwinian theories. This result corroborates what has been found by Ferraz-Mello et 
al. (2009): The second-degree expansion of the components of the force obtained with the lag time 
introduced by Mignard (1979) is equivalent to that obtained with Darwin standard theory when lags are 
assumed proportional to frequencies. The last equation can be written in terms of Mignard's constant 
time lag r via 5 — nr. 



5 Singer's modification 

Singer (1968) also proposed a modification in MacDonald's model to take into account the nonunifor- 
mity of /. He used MacDonald's equations, but considered that the distance from the vertex of the 
Earth tide to the Moon is variable, and used 

«-r(fl-|). (9) 

Since r = const, Singer's extension is sometimes considered as a pionneer example of theory with con- 
stant time-lag; however, when the patch proposed in eqn. Q is introduced in MacDonald's equations, 
the result does not correspond to a constant time-lag rheology as do, for instance, the Williams- 
Efroimsky equations. As in the original version of MacDonald, no law exists ruling the formation of 
the lags in Singer's tide components. The correction of the geometric lag nonuniformity entails, cor- 
rectly, the argument 2(1 — 0) + 2(n — Q)t in the semi-diurnal terms of Singer's U. However, in the other 
tide components the arguments do not follow any law relating lags and frequencies. For instance, we 
mention the monthly tide arguments (3^ — 20 + 3nr — 3Qt) and (£ — 20 + 2nr — J7r), the argument 20 
(without any lag), etc. Had Singer's patch been enough to endow MacDonald's model with Darwin's 
rheology, all terms would have the form k£ — k'4> + (kn — k! Q)t. 
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6 Conclusions 

The figure of the tidally deformed Earth in MacDonald's theory has its symmetry axis at a constant 
distance from the Moon. In that theory, the fact that the Moon is on an elliptic orbit only affects 
the ratio of the figure axes. The only variable is the prolateness of the spheroid, which is inversely 
proportional to the cube of the Earth-Moon distance. As a consequence, MacDonald's equations do 
not allow to classify his theory according with a rheology. The lags of the harmonic tide components 
only follow the geometrical constraints of this spheroid. They do not obey any law and are unrelated 
to the frequencies of the components. The patch proposed by Singer (1968) corrects the arguments of 
the terms describing the semi-diurnal tide, but for the remaining terms any law is followed. In order 
to get one constant lag-time theory, it is not sufficient to correct the nonuniformity of S, but it is 
also necessary to introduce Mignard's operational hypothesis (Mignard, 1979) after which the actual 
potential on one point on the Earth, at a given time, is equal to the potential which would be acting 
on it at a time r before the given time, would the Earth be inviscid. This operational hypothesis is 
also adopted in the modified MacDonald theory presented by Williams and Efroimsky (2012) and is 
instrumental to endows it with Darwin's rheology. 
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